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1 Introduction 

Nowadays, the dynamical description of strongly interacting systems out of equilibrium is 
dominantly based on transport theories and efficient numerical recipies have been set up 
for the solution of the coupled channel transport equations |]I|, ^, ^, ^, |^, |^, 0| (and Refs. 
therein). These transport approaches have been derived either from the Kadanoff-Baym 
equations f^ in Refs. [P, |l^, [ll], |12|, |l^] or from the hierarchy of connected equal-time 



Green functions M, IT5[ in Refs. H, ITB, IT^ by applying a Wigner-transformation and 



restricting to first order in the derivatives of the phase-space variables (X, P). 



However, as recognized early in these derivations [11, 16], the on-shell quasiparticle 
limit, that invoked additionally a reduction of the 8iV-dimensional phase-space to 7A^ 
independent degrees of freedom, where N denotes the number of particles in the system, 
should not be adequate for particles with high collision rates (cf. also Refs. [0, p!9||). 
Therefore, transport formulations for quasiparticles with dynamical spectral functions 



have been presented in the past [|T2, 12 1 providing a formal basis for an extension of 



the presently applied transport models denoted by BUU/VUU g, |, ||, 0, |22|, QMD 



|3|, ^ or its relativistic versions RBUU g], UrQMD §, ART 0, ARC O or 
HSD0,|1. 

Recently, the authors have developed a semiclassical transport approach on the basis 
of the Kadanoff-Baym equations that includes the propagation of hadrons with dynamical 
spectral functions [^. This approach has been examined for nucleus-nucleus collisions 
at GANIL energies and the off-shell propagation of nucleons and A's has lead to an 
enhancement of the high energy proton spectra as well as to an enhancement of high 
energy 7-rays rather well in line with related experimental studies. The question now 
arises whether such phenomena will also be encountered at SIS energies or AGS energies, 
i.e. at laboratory energies that are higher by one or two orders of magnitude, respectively. 

The paper is organized as follows: In Section 2 we will briefly review the general- 
ized transport equations on the basis of the Kadanoff-Baym equations [§, extend the 
test-particle representation to the general case of momentum-dependent self energies and 
specify the collision terms for bosons and fermions. Furthermore, we will present the mass 
differential cross sections for mesons from NN and vrA^ collisions in the medium as well as 
the approximations for their 'coUisional broadening'. The results for particle production 
in nucleus-nucleus collisions from the off-shell approach will be discussed in Section 3 in 
comparison to the quasi-particle limit as well as to the experimental data available. A 
summary and discussion of open problems concludes this study in Section 4. 

2 Extended semiclassical transport equations 

In this Section we briefly recall the basic equations for Green functions and particle self 
energies as well as their symmetry properties that will be exploited in the derivation of 
transport equations in the semiclassical limit. 



The general starting point for the derivation of a transport equation for particles with 
finite width are the Dyson-Schwinger equations for the retarded and advanced Green func- 
tions 5'''^*, S"''^'" and for the non-ordered Green functions S*^ and 5^ pOl- I^ the case of 



scalar bosons - which is considered in the following for simplicity - these Green functions 
are defined by 
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They depend on the space-time coordinates x, y as indicated by the indices -xy The Green 
functions are determined via Dyson-Schwinger equations by the retarded and advanced 
self energies S^^*, H"-'^'" and the collisional self energy E*^: 
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Equation (H) is the well-known Kadanoff-Baym equation. Here Sqx denotes the (negative) 
Klein-Gordon differential operator which is given for bosonic field quanta of (bare) mass 
Mq by Sqx = —{d!^d^ + Mq), Sxy represents the four-dimensional (5-distribution Sxy = 
S^^\x — y) and the symbol indicates an integration (from — cxd to cxd) over all common 
intermediate variables (cf. |l30| ). 

For the derivation of a semiclassical transport equation one now changes from a pure 
space-time formulation into the Wigner-representation. The theory is then formulated 
in terms of the center-of-mass variable X = [x + y)/2 and the momentum P, which is 
introduced by Fourier-transformation with respect to the relative space-time coordinate 
{x — y). In any semiclassical transport theory it is furthermore assumed that the depen- 
dence on the mean space-time coordinates X of all functions is rather weak. Therefore, 
in the Wigner-transformed expressions only contributions up to the first order in the 
space-time gradients are kept. After carrying-out these two steps the Dyson-Schwinger 
equations (HH) become 
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The operator O is defined as [|T2|, ^ 

1 / dFi dF2 dFi dF2 



{Fi}{F2} 



2 I ^X^' dP„ 



dP^ (9X/^^ 



which is a four-dimensional generalization of the well-known Poisson-bracket. Starting 
from (^ and (|^) one obtains algebraic relations between the real and the imaginary part 
of the retarded Green functions. On the other hand eq. (0) leads to a 'transport equation' 
for the Green function S^. 

We briefly recall the necessary steps: we separate all retarded and advanced quantities 
- Green functions and self energies - into real and imaginary parts, 
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The imaginary part of the retarded propagator is given (up to a factor 2) by the normalized 
spectral function 
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(10) 



while the imaginary part of the self energy corresponds to half the width Txp- By 
separating the complex equations (^ and (P) into their real and imaginary contributions 
we obtain an algebraic equation between the real and the imaginary part of S**"^*, 
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In addition we gain an algebraic solution for the spectral function (in first order gradient 
expansion) as 
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while the real part of the retarded propagator is given by 
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Furthermore, the (Wigner-transfornied) Kadanoff-Baym equation (|^ allows for the con- 
struction of a transport equation for the Green function S^. When separating the real and 
the imaginary contribution of this equation we find i) a generalized transport equation, 

0{p2 _ M^ - ReE^^'p}{S<p} ~ 0{J:<p}{ReSrp} 

and ii) a generalized mass-shell constraint 

[P' - Ml - i?eS-^]^<p - T.<pReSfp 

= lo{^lp}{Axp} -lo{rxp}{S^p}. (15) 

We note that so far no reduction to the quasiparticle limit has been introduced. Be- 
sides the drift term (i.e. OjP^ _ Ml}{S<} = -P>'d^S<) and the Vlasov term (i.e. 
0{PeS'-^*}{^<}) a third contribution appears on the l.h.s. of (0) (i.e. 0{S<}{i?e5'^^*}), 
which vanishes in the quasiparticle limit and incorporates - as shown in |^ - the off-shell 
behaviour in the particle propagation. It is worth to point out that this particular term 
has been discarded in Ref. |^ and an off-shell propagation had to be introduced by 
hand via some mass- dependent real potential. To our best knowledge, such an auxiliary 
potential cannot be extracted from the Kadanoff-Baym equation (^). The r.h.s. of (|14|) 
consists of a collision term with its characteristic gain (~ S^S*^) and loss (~ S^S^) 
structure, where scattering processes of particles into and out off a given phase-space cell 
are described. 

Within the specific term (0{S^}{i?eS"''^*}) a further modification is necessary. Ac- 
cording to Botermans and Malfiiet |jTO| the coUisional self energy S*^ should be replaced 
by S*^ • r/A to gain a consistent first order gradient expansion scheme. The replacement 
is allowed since the difference between these two expressions can be shown to be of first 
order in the space-time gradients itself. Therefore it has to be discarded when appearing 
inside an additional Poisson-bracket. Furthermore, this substitution is required to get 
rid of the inequivalence between the general transport equation and the general mass 
shell constraint which is present in (|T^) and (|T^) [|TS[ but vanishes by invoking S^ ■ T/A. 



Finally, the general transport equation (in first order gradient expansion) reads |3( 
Axp Fxp 

= i [T,j^p Sxp — T^xp ^xp\ ■ (16) 



0{P'-M',-ReJ:%}{S<p} - ^0{Txp}{iP'-M',-ReJ:Tp)S^p} 

J- XP 



Its formal structure is fixed by the approximations applied. Note, however, that the 



dynamics is fully determined by the different self energies, i.e. ReT,^xp, Txp, ^xp ^^^ 
S^p that have to be specified for the physical systems of interest. 

2.1 Testparticle representation 

In order to obtain an approximate solution to the transport equation (|T6|) we use a 
testparticle ansatz for the Green function S^, more specifically for the real and positive 
semidefinite quantity 

N 

FxP =AxpNxP= iS<p ~ Y. ^^'\^ - Mt)) S^'HP - m) S{Po-e,{t)).{17) 

i=l 



Whereas so far we have briefly repeated the derivation from [^, we now extend the 
testparticle description to explicitly four-momentum-dependent self energies. The latter 
is not essential for the description of low and intermediate energy heavy-ion reactions, 
but becomes necessary for relativistic dynamics. In the most general case (where the 
self energies depend on four-momentum P, time t and the spatial coordinates X) the 
equations of motion for the testparticles read 
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where the notation F(j) implies that the function is taken at the coordinates of the test- 
particle, i.e. F(i) = F(t,Xi(t),Pi{t),ei(t)). 



In ([T8| - P0D a common multiplication factor (1 — C(j)) appears, which contains the 
energy derivatives of the retarded self energy 
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It yields a shift of the system time t to the 'eigentime' of particle i defined by tj = t/(l — 
C(j)). As the reader immediately verifies, the derivatives with respect to the 'eigentime', 
i.e. dXi/dti, dPi/dti and dei/dti then emerge without this renormalization factor for each 
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testparticle i when neglecting higher order time derivatives in hne with the seniiclassical 
approximation scheme. The role and the importance of this correction factors will be 



studied in detail for a four-momentum-dependent 'trial' potential in Section p74 We 
note that for momentum-independent self energies we regain the transport equations 
as derived in [Q; only in case of particles with a vanishing vacuum width Ty = and 
Txp ~ Pb (baryon density) these equations reduce to the 'ad hoc' assumptions introduced 
in Ref. [^. Furthermore, in the limiting case of particles with vanishing gradients of the 
width Txp these equations of motion reduce to the well-known transport equations of the 
quasiparticle picture. 



Following Ref. |3^ we take M^ = P^ — ReY7'^^ as an independent variable instead of 



Pq; which then fixes the energy (for given P and M^) to 

^ret 
■'XPM'2 



P^ = P^ + M^ + ReJ:fp^, . (22) 



Eq. ( pOD then turns to 

dMl ^ Ml - Ml ^ 
dt r(j) dt ^ ' 

for the time evolution of the test-particle i in the invariant mass squared as derived in 
Ref. ig. 



We briefly comment that in Ref. [^ we have added a term ~ dV{i)/dt in the equation 
of motion for the test-particle momenta in order to achieve strict energy conservation 
for each test-particle; however, within the subsequent approximations introduced in the 
actual transport calculations for energetic nucleus-nucleus collisions such a term did not 
lead to noticable effects within the numerical accuracy achieved. We thus discard such a 
term in the present formulation and investigation. 

2.2 Model simulations for the momentum-dependent transport 
equations in testparticle representation 



To demonstrate the physical content of the equations of motion for testparticles ([l8| - pOD we 
perform an exploratory study with a momentum-dependent trial potential. The potential 
is chosen of the type: 



.ret V_ V{P,,P) J WiPo,P) 



2 l + exp{(|rl-P)/ao} \^1 + exp{(|f| - P)/ao} 2) ^ ^ 

with a constant (but finite) vacuum width Ty. While the spatial extension of the potential 
is given (as in ||30|) by a Woods-Saxon shape (with parameters P = 5 f m and Oq = 0.6 
fm) its momentum dependence for the real as well as for the imaginary part is introduced 

by 

F(Po,P) =Cv ^ ^, W{Po,P) =Cw ^ ^. (25) 

K-{Pi-P^) Kr-{Pi-P') 



Here the constants Cy iCw) give the 'strength' of the complex potential while Ky {^w) 
play the role of cutoff-parameters. Due to the structure in the denominator of ( P5| ) the 
momentum-dependent part of this potential is explicitly Lorentz-covariant. 

In our simulation we propagate the testparticles with different initial mass parameters 
Mj, which are shifted relative to each other by rv'/(20 GeV) around a mean mass of 1.0 
GeV. To each testparticle a momentum in positive ^-direction is attributed so that all 
of them have initially the same energy Pq = 2.0 GeV. All testparticles are initialized on 
the negative 2;-axis with {\Xi{t = 0)| ~ 15 fm) and then evolved in time according to the 
equations of motion (p!^-pO|). 

In our first simulation we consider a purely imaginary potential with a strength of 
Cw = 0.6 GeV^ and a cutoff-parameter A^y = 2.0 GeV. The evolution in energy Pqj; 
momentum P^i and in the mass parameter M, for all testparticles is shown in Fig. 1 
(upper part) as a function of z(t). When the testparticles enter the potential region, their 
momenta and mass parameters are modified. As already shown in [^ the imaginary 
potential leads to a spreading of the trajectories in the mass parameter M, which in turn 
reflects a broadening of the spectral function. The relation between the imaginary self 
energy and the spreading in mass is fully determined by relation (^). Since we have 
chosen a potential with no explicit time dependence the energy of each testparticle is a 



constant of time (cf. Ref. |^). According to the explicit momentum dependence of 
our 'trial' potential each single testparticle is affected with different strength. Since the 
imaginary potential is strongest for small momenta (which correspond to the highest lines 
in the lower graph of Fig. 1) the momentum and mass coordinates of those testparticles 
are changed predominantly that are initialized with the lowest momenta (i.e. with the 
largest masses). As a result one observes a rather asymmetric distribution in the mass 
parameters (and in the momenta) in the potential zone. This is different from the studies 
of 1^ where the investigated momentum-independent potential yields a nearly equidistant 
spreading of the mass trajectories. For z{t) 3> R the mass and momentum coordinates of 
the testparticles return to the proper asymptotic value. 

In the second example we allow for an addititonal real part of the self energy. The 
calculation is performed with the parameters Cy = —0.3 GeV^, Cw = 0.6 GeV^ and 
Ay = Aw = 2.0 GeV. The momentum-dependent real part (lower part of Fig. 2) causes 



as also observed in ||3^ - an additional shift of the testparticle momenta. Since the 



real part of the potential is larger for small initial momenta, these testparticle momenta 
are shifted up somewhat more than for particles with larger momenta. This gives rise 
to a reduction of the asymmetry which was introduced by the momentum-dependent 



imaginary part of the self energy (upper part of Fig. 2). As in |30] the mass parameters 
of the testparticles are only weakly influenced by the real part of the potential. 

We, furthermore, study the implications of the correction term (1 — C(i))~^ using the 
same imaginary potential as in Fig. 1. The time evolution of the correction factor for each 
testparticle i is displayed in the lower part of Fig 3. While it is > 1 for large initial mass 
parameters Mj > Mq, it is < 1 for mass parameters Mj < Mq. In the upper part of Fig. 



3 the momenta of the testparticles are shown for two calculational hmits: in the first one 
the correction term is taken into account (as in the previous calculations), while in the 
second one the corrections due to the energy dependence of the retarded self energy are 
neglected. However, the calculations with and without the correction factor exhibit only 
a very small difference in the testparticle momenta, which even cannot be distinguished 
within the resolution of Fig. 3. The same holds for the mass parameters Mi which are not 
displayed here since they provide no new information due to energy conservation. We thus 
conclude that the particle trajectory is not very sensitive to these correction factors, since 
the correction term - when displayed as Pz{z) in phase-space - leads only to a rescaling 
of the 'eigentime' of the testparticles as pointed out before. 

2.3 Collision terms 

The collision term of the Kadanoff-Baym equation can only be worked out in more detail 
by giving explicit approximations for S^ and S^. A corresponding collision term can 



be formulated in full analogy to Refs. [y, |T6|, e.g. from Dirac-Brueckner theory, and 
implementing detailed balance as 

hoii{X, P, M^) = Tr2TrsTnA{X, P, M^)A{X, A, M^)A{X, Ps, M|)A(X, P^, m^) 
|r((P,M2) + {P2,M^) ^ {P3,M',) + {P,,Ml))\Xs S'^'\P + P2 - Ps - P,) 

[^XPaAq^XPiM^fxPlVpfxPiMl ~ ^XpM^ ^XPiM^ f'xPsMl fxP^M'il (26) 

with 

and 1] = ±1 for bosons/fermions, respectively. The indices A,S stand for the antisym- 
metric/symmetric matrix element of the in-medium scattering amplitude T in case of 
fermions/bosons. In eq. (^) the trace over particles 2,3,4 reads explicitly for fermions 

Tr2 = Y. TTT^T-J dP2 1^ ' , (28) 



-^2(27rr-^ 2VP| + M; 



'2 



where cr2, T2 denote the spin and isospin of particle 2. In case of bosons we have 

^--Epi)./<^'P#, (2.) 

since here the spectral function As is normalized as 

J^Ab{X,P) = 1 (30) 



whereas for fermions we have 
dPo 



J^A,{X,P) = 1. (31) 



We mention that the spectral function Ap in case of fermions in ( PBD is obtained by 
considering only particles of positive energy and assuming the spectral function to be 
identical for spin 'up' and 'down' states. In general, the spectral function for fermions 
AajsiX, P) is a Dirac-tensor with ajS denoting the Dirac indices. It is normalized as 

rIP 

-^A^,{X,P) = {^Xp, (32) 

which implies 

\T.l^i^'Mx,P)U = l. (33) 

Now expanding A in terms of free spinors Us{P) (s=l,2) and Vs{P) as 

(i),^ = Y. ^s{PM)p Al, Ur{P, M), + Vs{PM)f^ ^r? ^P, M), (34) 

r,s=l 

one can separate particles and antiparticles. By neglection of the antiparticle contributions 
(i.e. A'^'P = 0) and within the assumption that the spectral function for the particles is 
diagonal in spin-space (i.e. A^^ = SrsA^) as well as spin symmetric, one can define Ap as 

Ap^A{ = Al. (35) 

Neglecting the 'gain-term' in eq. (p6|) one recognizes that the coUisional width of the 
particle in the rest frame is given by 

T,ai{X, P, M') = Tr^TrsTr, |T((P, m') + (P^, M^) ^ (P3, M|) + (P^, Ml))\Xs (36) 
AiX, P2, Mi)AiX, P3, Mi)AiX, P4, M!) 5\P + P2 - P3 - P4) N^p^^,, f^p,^,,, f^p^^,, , 

where as in eq. (|26| ) local on-shell scattering processes are assumed for the transitions 
P + P2 — > P3 + P4. We note that the extension of eq. (|26|) to inelastic scattering 
processes (e.g. A^A^ -^ A^A) or [ttN —>■ A etc.) is straightforward when exchanging 
the elastic transition amplitude T by the corresponding inelastic one and taking care of 
Pauli-blocking or Bose-enhancement for the particles in the final state. We note that for 
bosons we will neglect a Bose-enhancement factor throughout this work since their actual 
phase-space density is small for the systems of interest. 

For particles of infinite life time in vacuum - such as protons - the coUisional width 



(|36D has to be identified with twice the imaginary part of the self energy. Thus the 
transport approach determines the particle spectral function dynamically via ( p^ ) for 
all hadrons if the in-medium transition amplitudes T are known in their full off-shell 
dependence. Since this information is not available for configurations of hot and dense 
matter, which is the major subject of future development, a couple of assumptions and 
numerical approximation schemes have to be invoked in actual applications. 
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2.4 Numerical realization 

As in Ref. |^D| the following dynamical calculations are based on the conventional HSD 
transport approach |^, ^ - in which ReHxp is specified for the hadrons - however, the 



equations of motion for the testparticles are extended to (08,19,23). Whereas for energies 



up to 100 A MeV (GANIL energies) essentially the nucleon degrees of freedom were 



important |]30| we have to specify the actual 'recipies' for nucleon and meson resonances 
involved in the following calculations. 

As a first approximation we only consider reactions with binary channels which for 
convenience can be described in their center-of-mass (cms) frame. The collisions of nu- 
cleons (as well as all hadrons) are described by the closest distance criterion of Kodama 
et al. [^ : a collision of two particles takes place only if their distance in the individual 
cms is small enough, i.e. 



|Xi-X2|< V^dPi-PaD/vr, (37) 

where a denotes the total cross section of the process, which is written here as a function 
of the three-momentum difference in the cms. In the case of nucleon-nucleon collisions 
the Cugnon parametrization ||33[ for the in-medium NN cross section da / dVL[\fs') is used 
by identifying (in the A^A^ c.m.s.) 

s' - 4m% = s- 4M^ = 4P^ (38) 

where mi\[ is the nucleon vacuum mass, M the actual off-shell mass and vs' the invariant 
energy of a nucleon-nucleon collision in the vacuum with the same cms-momentum P. The 
final nucleon states are selected by Monte-Carlo according to the local spectral function 
determined by the collisional width TcouiX, P), while the angular distribution in the cms is 
taken the same as for on-shell nucleons. This recipe for off-shell nucleon-nucleon scattering 
is practically an ad-hoc assumption and has to be controlled by off-shell matrix elements 
of the nucleon Brueckner T-matrix in the medium. In Section 3 we will also investigate 
an alternative recipe to demonstrate the model dependence of our actual results. 

If two particles have approached sufficiently close, furthermore, the time of the collision 
in the calculational frame has to be specified. Within a relativistic treatment this is a 
nontrivial task since 'simultaneity' is a non-invariant property. Therefore, it is necessary 
to take into account the time coordinates ti and ^2 for both participants of the collision 
separately. By determining the collision time (as the time of the closest approach) in 
the individual cms of each colliding particle, a causality respecting time ordering of all 
collisions in the calculational frame can be achieved. Without representing the explicit 
formulae we refer the reader for a detailed description to eqs. (8)-(ll) of Kodama et al. 
p2| which are implemented in the HSD transport approach. 

In order to determine the cross section a in { ^7\ ) let us consider a binary reaction of 
hadrons characterized by Pi,Mf and P2,^2 suppressing all internal quantum numbers. 
In case of a single final state, i.e. a resonance R for meson-nucleon or meson-meson 
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scattering, energy and momentum conservation fixes the four-momentum of the final state 
completely. The cross section for the production of the resonance R in the center-of-mass 
system with invariant energy M then is given by the Breit-Wigner cross section 



a{M^ 



C. 



s,i 
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Pin 



M 



2p2 ID 



(M2 - M^f + M2r 



2 • 
tot 



(39) 



where pf^ is the squared three-momentum of a particle in the entrance channel, Bin the 
branching ratio to the resonance R and 



tot 



Ty{M'') + T,,u{M^ 



r 



XP 



2Pn 



(40) 



In (^OD Ty and Tcou denote the vacuum and collisional particle width, respectively; 
it is related to Vxp = -2JmE''^* by Txp = 2Po Ttoi H- In (|9|) the factor Cs,i 
is the usual spin/isospin factor determined by the spin/isopin in the entrance chan- 
nel and the resonance properties. The vacuum decay width for the A is taken in the 
parametrization of p2| corrected by the ratio of the two-body phase-space integrals 
R2{^/s,Ml,M2)/R2{^/s,M'^,M2). Here the two-body phase-space integral is given by 
Mi 



R2( 3,171^,171^'^ 



^l5 



IT 



X'^[s,m", m\ 



^l 



with 



A(s, ml, ml 



s — {mi + 1712^ s — (mi — m.2) 



(41) 



(42) 



For the iV(1440) we assume the same parametrization of the width as in Ref. [^ corrected 
again by the ratio of two-body phase-space integrals. The same strategy is taken for the 
A^(1535) [^] which is important for rj production, rescattering and absorption. Higher 
baryon resonances (as in Refs. [^, 



36[) are not considered in this study since they are 
not seen experimentally in the photon absorption experiments in Frascati even on light 
nuclei |]3^ . 

In case of binary exit channels the final masses M3,M4 are selected by Monte-Carlo 
according to the local spectral functions with width Ttot- The A^A^ -^ NR transitions are 
simulated by i) correcting the vacuum width according to the available phase-space in the 
final channel and ii) by adding the collisional width Tcoii- These assumptions presently 
are also hard to control by fully microscopic calculations and serve as a guide for the 
effects to be investigated below. 

In case of meson production by off-shell baryon-baryon or meson-baryon collisions we 
either have 2 (e.g. nN -^ K+A/S), 3 (e.g. A^A^ -^ K+AN or K+EN) or 4 particles (e.g. 
A^A^ — >■ NNK^K~) in the final channel. Since the final mesons may be off-shell as well, 
one has to specify the corresponding mass-differential cross sections that depend on the 
entrance channel and especially on the availabe energy ^/s in the entrance channel. 
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We start with the exphcit parametrizations for meson {m) production cross sections 
given in Ref. ^ for on-shell mesons as a function of the invariant energy y/s in case 
of nucleon-nucleon or pion-nucleon coUisions, i.e. <7NN^NNm{y/s) or cr-^N-^mNiy/s), re- 
spectively, that are well controlled by experimental data. Far above the corresponding 
thresholds the mass differential cross sections are approximated by 

^^^ '^ = aNN^mNN{\Js - ^ sl) A^{M %T tot) , (43) 

where Am^M"^ ,T tot) denotes the meson spectral function for given total width Ttot that 
is normalized to unity by integration over dM'^. In (|43| ) the threshold energy -^/sj = 
Mo + Ml + M2 depends on the masses of the hadrons in the final channel, i.e. Mq, Ml 
and Mg. Actual events then are selected by Monte-Carlo according to (|43|). Close to 
threshold v^, i.e. for ^/s - Mq - Ml - M^ < 2Ttot, where Ml,M^ denote the final 
off-shell masses of two nucleons, Mq the meson pole mass and Ttot its total width, the 
differential production cross section is approximated by a constant matrix element squared 
|Mmp times available phase-space, 

da^NjMV^s) ^ |^^|2 ^^(^2^ r,„,) R,{s, M\ Ml\ Mt). (44) 

The matrix element \Mm\ then is fitted to the on-shell cross section close to threshold. 
In (^^ the function R^ denotes the 3-body phase-space integral in case of a mNN final 
state and is given by pi 



2 2 2\ /"(v^-™!)^ dS2 ,1/ 2\ a!/ 2 2\ 

R3{s,m^,m2,m^) = / A2(s2, s,mj Aa (52,^2, mg) . (45) 

J(mo+m-i)'^ So 



The same recipe is used for binary mN channels by replacing R^ with the 2-body phase- 
space integral R2. 

In case of 4 particles in the final state, e.g. in the channel NN -^ KKNIN2, where 
the A^*'s denote off-shell nucleons, the differential cross section is approximated by 

d^Pid^P2 d^ps d^p4 

0"BB^Ar7VMiAf2(v^)T7^"^^V^^(^l + ^2 - Pi - P2 - Ps - Pi), 

where R4 denotes the 4-body phase-space integral [Q. Similar strategies have been 



exploited in case of subthreshold pp production in proton-nucleus and nucleus-nucleus 



collisions in Refs. [^, [39|] 



The resulting cross sections for K production from A^A^ and n'^p collisions are dis- 
played in Fig. 4 as a function of y/s for different collisional width Tcou= 0, 50, 100, 150, 
200 MeV while keeping ^ = Mk + Mn + M^ or ^ = Mk + Ma, respectively. With 
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increasing width Tcou the subthreshold production of mesons becomes enhanced consid- 
erably relative to the respective vacuum cross section, but the absolute magnitude stays 
small below threshold even for Tcoii = 200 MeV. 

As a next step we have to fix the coUisional width of the hadrons in the nuclear medium 
which enters the spectral function A{X,P,M'^) as well as the differential cross sections 
(^31) ■ According to ( |36D the coUisional width is explicitly momentum- (and energy-) 
dependent. Whereas in Ref. [^ we have employed momentum-independent coUisional 



broadening, this is no longer adequate for relativistic systems. We thus evaluate Tcou for 
each particle in a finite cell in coordinate space according to (p6D, however, discard the 
explicit dependence on the energy Pq. This approximation implies that the correction 
factors (1 — C(j))~^ in the testparticle equations of motion (18) - (20) are equal to 1 and 
all particles can be propagated with the same system time. 

In case of kaons, antikaons or p mesons at SIS energies we treat the latter perturba- 



tively as in Ref. |^, i.e. each testparticle achieves a weight Wi defined by the ratio of 
the individual production cross section to the total ttB or BB cross section at the same 
invariant energy. Their propagation and interactions are evaluated as for baryons and 
pions, however, the baryons (pions) are not changed in their final state when interacting 
with a 'perturbative' particle. The actual coUisional width then is approximated by 

ri ~ Sj Vij(^ij . ,„s 

■■- coll ~ 7i y^ ; l^ ' J 

where the sum over j runs over all baryons in the local cell, Vij is the relative velocity in 
the meson-baryon cms, 7j is the Lorentz-factor of the particle with respect to the local 
rest frame of the baryons and aij their total cross section at invariant energy ^/s--. Note, 
that in (|47| ) the final state Pauli-blo eking has been neglected for the baryon which should 
be reasonable at the high bombarding energies of interest here. 

Apart from the description of particle propagation and rescattering the results of 
the transport approach also depend on the initial conditions, Xj(0), Pj(0), Mf (0). In 
view of nucleus-nucleus collisions, i.e. two nuclei impinging towards each other with 
a laboratory momentum per particle Piab/A, the nuclei can be considered as in their 
respective groundstate, which in the semiclassical limit is given by the local Thomas- 
Fermi distribution |^ . Additionally the virtual mass Mf for nucleons has been determined 
by Monte-Carlo according to the Breit-Wigner distribution (|12D assuming an in-medium 
width To = 1 MeV. We mention that varying Fq from 1-5 MeV does not change the 
results to be presented in Section 3 within the statistics achieved. For the vacuum width 
of stable hadrons we have used Ty = 1 MeV which implies that nucleons propagating to 
the continuum in the final state of the reaction achieve their vacuum mass on the 0.1 % 
level. 
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3 Nucleus-nucleus collisions 

Our concrete applications we first carry out for nuclear reactions at SIS energies (1-2 
A GeV) that have been analysed within conventional transport models to a large extent 
(cf. Ref. [0] and Refs. cited therein). 

In view of Eq. (p3D we present for some randomly chosen testparticles i the off-mass- 
shell behaviour Mf{t) — Mq as a function of time in a central collision (6 = 1.5 fm) in Fig. 
5 for Au + Au at 1 A GeV. It is seen that during the collision of the nuclei from t ^ 7 - 25 
fm/c the off-shellness of baryons reaches up to 0.8 GeV^, however, the nucleons become 
practically on-shell for t > 35 fm/c. The individual sudden high mass excitations and 
subsequent decays correspond essentially to A and N(1440) baryons. Nucleons in their 
decay may be off-shell, but propagate again to their on-shell mass in the continuum. 

The baryon spectral function is shown in Fig. 6 as a function of the invariant mass for 
the latter reaction at times of 0, 5, 10, 20, 40 and 60 fm/c. Apart from a broadening of 
the nucleon spectral function at the initial time steps one observes that the high mass tail 
is completely covered by the A and N(1440) excitations. This is different from the results 
in |]3y] at GANIL energies since in the latter case the A excitation was dynamically 



suppressed and the high mass tail of the spectral function dominated by nucleons. In 
fact, in 1 A GeV Au + Au collisions the resonance high mass spectrum in the off-shell 
calculations is only slightly enhanced as compared to the on-shell calculations (without 
explicit representation). Note that at t = 60 fm/c all resonances have decayed and the 
nucleons have become on-shell again. 

The dominance of the resonances for the high mass tail can also be seen in the dif- 
ferential collision number dNcoll/d^/s for baryon-baryon collisions which is displayed in 
Fig. 7. Here the dashed histogram stands for the on-shell result while the solid histogram 
gives the off-shell distribution that extends well below the two-nucleon threshold. The 
high energy tail in this distribution dominantly arises from A^A and AA or A^A^(1440) 
reactions, which are similar in both calculations. We thus find only a small enhancement 
in the high energy tail for the off-shell case relative to the on-shell limit. 

The latter observation can also be made in the transverse proton momentum spectra 
1/pTdNp/dpT (Fig. 8) for Au + Au at 1 A GeV (6=1.5 fm) where the on-shell propagation 
(dashed histogram) practically leads to the same result as the off-shell propagation (solid 
histogram) except for the high momentum tail. Since the elastic A^A^ cross section is 
taken here as a function of the momentum difference in the cms (cf. Section 2.4) one 
might worry if alternative prescriptions for aNN could change the results. In this context 
we have performed calculations using the Cugnon parametrization <7]\fj\f{y/s) with ^/s 
denoting the actual invariant energy of the off-shell nucleons and adopting a^N = 55 mb 
for invariant energies below 2 rriN, where ttln denotes the nucleon vacuum mass. The 
transverse momentum spectra according to this recipe (for the same reaction) are shown 
in Fig. 8 in terms of the dot-dashed histogram, which coincides with the 'default' off-shell 
result (solid histogram) within the statistics. 
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3.1 Meson production at SIS energies 

In principle also the pions should be propagated with a dynamical spectral function since 
their coupling to nucleons is very strong. In fact, pion collision rates in these reactions 
lead to a coUisional width which is much larger than the pion mass itself. A straight 
forward selection of pion masses (e.g. in the A decay) according to such a dynamical 
spectral function leads to M^ < which implies that such testparticles become acausal, 
i.e. move with velocities /3 > 1. In order to avoid such inconsistencies we treat pions 
on-shell throughout this study. We note, however, that microcausality is an essential 
issue that should also survive in transport approximations. This is practically done by 
the explicit requirement Mf > in the transport calculation, but not yet inherent in eqs. 

(Hl-il- 



Since the authors of Ref. ||31| claim a large enhancement in the pion yield for Au + Au 
at 1 A GeV and especially in the high energy tails of the pion spectrum when describing 
nucleon off-shell propagation in their Monte-Carlo simulation, we show in Fig. 9 (l.h.s.) 
the inclusive differential vr"*" spectrum for this system at diab = 44 ± 4° within our off-shell 
approach (solid histogram) in comparison to the on-shell limit (dashed histogram) and 
the experimental data from the KaoS collaboration |^0|]. As in Ref. |^ the vr"*" spectrum 
is described quite well within the on-shell HSD approach. As also seen from Fig. 9 there 
is only a very slight enhancement of the high momentum tail in the pion spectrum for 
our off-shell calculation which is still compatible with the experimental spectrum from 



Ref. [|y] within the statistical errors. Note that the fluctuations in the histograms with 
respect to an average exponential spectrum provide some information about the statistical 
accuracy. Thus our off-shell approach - based on the Kadanoff-Baym equation (|16|) - is 
not in conflict with the experimental data contrary to the model from Ref. ||3T[| . 

Since kaons couple only weakly to nucleons and are not absorbed at low energies 
their coUisional width is rather small such that they may be treated on-shell to a good 
approximation. In their differential production cross section we thus only can test the 
effects from the off-shell propagation of baryons. The inclusive K~^ spectra at 6iab = 44±4° 
for Au + Au at 1 A GeV are shown in Fig. 9 (r.h.s.) where we compare our off-shell results 
(solid histogram) with the on-shell limit (dashed histogram) and the experimental data 
from Ref. |^ (full circles) and Refs. |^2|, |^ (open circles) where the latter differ on 
average by a factor of 2. The off-shell calculations give a slightly higher K~^ yield than the 
on-shell calculations, however, are still well within the error bars of the presently available 
data. We mention that we have performed the calculations without any kaon potential; a 
slightly repulsive kaon potential will drop the K~^ yield and harden the spectrum slightly 
as discussed in Ref. P[ . 

The relative increase of the K~^ spectrum by about 50% in the off-shell calculation 
with respect to the on-shell result should also be compared with the relative sensitivity to 
the incompressibility K of the nuclear equation-of-state (EoS) [^. According to the early 



relativistic calculations by Lang et al. [|g] the inclusive K~^ yield in Au + Au collisions 
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at 1 A GeV is enhanced by about a factor of 2 when decreasing the incompressibihty K 
from 380 MeV to 200 MeV. When restricting to incompressibihties 200 MeV < iT < 300 
MeV this relative sensitivity reduces to a similar change in the K^ yield for Au + Au at 
1 A GeV as obtained from the off-shell dynamics relative to the on-shell treatment. It 
thus appears questionable if the incompressibihty of the EoS can be determined from the 
systematics of K^ spectra alone since also the K^ in-medium potential potential is not 
yet determined very accurately. 

We continue with pion production in Ni + Ni collisions at 1.8 A GeV since for this 



system also K^ and K~ spectra have been measured ||2, 50, J3|. The inclusive differential 
cross section for tt"'" from Ni + Ni collisions at 1.8 A GeV is shown in Fig. 10 as a function 
of the pion momentum in the cms in comparison to the data of the KaoS Collaboration 



taken at Qiab = 44 ± 4° [H2l H3] . The calculation with an off-shell propagation of baryons 



is represented by the solid histogram while the on-shell limit is displayed in terms of the 
dashed histogram. Within the numerical accuracy both results almost coincide which in 
view of Fig. 7 is not an exciting surprise. Both results, furthermore, are in a reasonable 
agreement with the measured spectrum (full squares) which suggests that no 'unphysical' 
approximations have been introduced within the 'recipies' for the off-shell transition cross 
sections. 

As might have been anticipated from the studies before, also the inclusive K~^ spectra 
from Ni + Ni at 1.8 A GeV - plotted as a function of the kaon momentum in the cms - 
is only weakly sensitive to the off-shell propagation of baryons as shown in Fig. 11 since 
the differential ^/s distribution for baryon-baryon collisions does not differ very much and 
also the pion-baryon production channel is similar (except for the high momentum tail). 
The latter fact one might also extract from the low sensitivity of the pion spectrum to an 
off-shell propagation in Fig. 10 for this reaction. We note that the kaons again have been 
propagated without any in-medium potential, which should be slightly repulsive in line 
with Refs. [^ |^, |5T|. As shown in Ref. [Q such a repulsive potential will suppress the 
kaon spectra especially at low momenta. Since we have employed the same production 
cross sections as in Ref. [Q we again reproduce the data from the KaoS Collaboration, 



taken at Qiab = 44 ± 4°, best without a kaon potential, where the off-shell calculation 
seems to be in even better agreement with the data. We mention that the 'theoretical' 
error bars in Fig. 11 are of statistical nature only and computed as ±0.5Sp/ JNp where Sp 
is the calculated spectral point in the actual momentum bin and Np denotes the number 
of events contributing to this momentum bin. Though the error bars between the on-shell 
(open triangles) and off-shell calculation (full triangles) almost overlap, we point out that 
the off-shell calculation gives a slightly harder spectrum. 

Note that the production channel A^A — >■ NK~^Y, where Y denotes a hyperon, as well 
as the AA -^ K^NY channel is not known experimentally and simple isospin factors 
as extracted from pion exchange |4^ might not be appropriate. Though there are some 



recent efforts to resolve this uncertainty within extended boson exchange models |55], the 



latter models will hardly be tested experimentally. This general uncertainty has to be 
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kept in mind when comparing transport calculations to experimental kaon spectra. 

We step on with the production of antikaons in the same reaction. Antikaons couple 
strongly to nucleons and thus achieve a large coUisional width in the nuclear medium. 
The calculations of Ref. |^2[|, which are based on a dispersion approach, give a coUisional 
width of about 100 MeV of K~ mesons at moderate momenta and nuclear matter density 
Po- Thus off-shell antikaons might be produced at far subthreshold energies (cf. Fig. 4), 
become asymptotically on-shell and thus enhance the K~ yield. Note, that this mechanism 
also might explain the K~ enhancement seen by the FRS, KaoS and FOPI collaborations 
|3|, |8], 3 in M + A^i reactions around 1.8 A GeV. 



Within our present 'modelling' of off-shell production processes, which are governed 
by phase-space close to threshold energies, we find an enhancement by about a factor ~ 2 
in the production of antikaons when treating baryons as well as antikaons off-shell. This is 
demonstrated in Fig. 12 where we show the K~ spectra as a function of their momentum in 
the cms for the off-shell (solid histogram, full triangles) and on-shell propagation (dashed 
histogram, open triangles). As in case of Fig. 11 the 'theoretical' error bars are statistical, 
only, and taken as iO.SS'p/w Ap where Sp and Np denote the actual value in the momentum 
bin and number of events, respectively. Both results underestimate the data from the FRS 
and KaoS Collaborations [|^, ^, ^ such that the conclusion in Refs. [|^, Q on attractive 



K self energies in the nuclear medium persists, though the off-shell calculations suggest 
somewhat smaller K~ potentials. We mention that the dispersion analysis of the K~ 



potential in nuclear matter in Ref. |^ also yields somewhat smaller antikaon potentials 
than that extracted in Refs. [^ ^ which is in line with our present finding. However, in 
order to obtain more model independent results on the antikaon self energy, precise data 
on antikaon flow from nucleus-nucleus collisions are urgently needed. 

3.2 AGS energies 

At AGS energies of ~ 11 A GeV the invariant energy s/s of the initial nucleon-nucleon 
collisions is about 5 GeV which implies that the nucleons are excited to continuum states 



denoted by strings - which decay according to phase-space [^ and fixed quark/diquark 



or s/{u^ d) ratios after a formation time Tp = 0.8 fm/c. Some details of the decay scheme 



are given in Ref. |5^. The FRITIOF 7.2 version - as implemented in the HSD transport 
approach - also includes the production of unstable particles of width F/j by Monte-Carlo. 
The width Th has been modified dynamically according to the actual total width at space- 
time point X which allows to simulate the production of hadrons with broad spectral 
functions from string decay in a straight forward manner. The low energy baryon-baryon 
and meson-baryon reactions are treated as for SIS energies (cf. Section 2.4). 

We here present only a single study for central Au + Au reactions at 11.3 A GeV {b < 
3.5 fm) and concentrate on the rapidity and transverse mass spectra of pions, kaons and 
antikaons. The results of our calculations are displayed in Figs. 13 and 14 for the off- 
shell calculation (solid histograms) and the on-shell calculations (dashed histograms) in 



comparison to the data from Refs. |]58| , |59| . Since at this bombarding energy the dynamics 
are dominated by continuum string excitations, a broadening of the spectral functions is 
found to play a minor role here. Note, that the K^ spectra for Au + Au collisions at 
AGS energies are underestimated in comparison to the data, which in Ref. [^ has been 
attributed to nonhadronic degrees of freedom or a partial restoration of chiral symmetry 
during the high density collision phase. 

4 Summary 

In this work we have employed the semiclassical off-shell transport approach from Ref. 



pO| , that in first order in the gradient expansion describes the virtual propagation of 
particles in the invariant mass squared M^ besides the conventional propagation in the 
mean-field potential (given by the real part of the retarded self energy), to analyse nucleus- 
nucleus collisions at SIS and AGS energies. Note that in conventional transport ap- 
proaches the imaginary part of the self energy is reformulated in terms of a collision 
integral and simulated by on-shell binary collisions, only. We here additionally account 
for the off-shell propagation of particles due to the imaginary part of the self energy 
in eqs. (|T8i l9,23) which describe the dynamical evolution of the particle spectral func- 



tion. The additional imaginary part of the self energy (or local collision rate Tcou) is 
determined by the collision integrals themselves and can be used in transport approaches 
without introducing any new assumptions or parameters provided that the off-shell tran- 
sition amplitudes T are known in the collision integral (^61) . In addition to Ref. Q] we 
have included momentum-dependent self energies for baryons and mesons in our present 
approach. However, a final answer to the role of off-shell hadrons in nucleus-nucleus col- 
lisions we still have to leave for future due to the 'ad hoc' assumptions involved in the 
treatment of the off-shell transition probabilities. 

We have presented dynamical calculations of the novel transport theory for nucleus- 
nucleus collisions at SIS and AGS energies where we can test its results in comparison 
to experimental data. We find that the off-shell propagation of nucleons practically does 
not change the rapidity distributions dN/dy and has only a minor effect on the trans- 
verse momentum spectra of protons within the statistics reached except for the very high 
momentum tails. The distribution of baryon-baryon collisions in the invariant energy 
\/i is found to be also enhanced only for high invariant energies since here the collisions 
with or between resonances - which are only slightly affected in their high mass spectrum 
- dominate the spectrum. Again except for high momentum tails we find no dramatic 
change in the pion and K^ spectra at SIS energies for Au + Au at 1.0 A GeV and Ni + Ni 
at 1.8 A GeV, our results being well in line with the data of the KaoS Collaboration. 
This no longer holds for the K~ spectra from Ni + Ni collisions at 1.8 A GeV which 
are enhanced by a factor of ~ 2 relative to the on-shell calculation within the statistics 
reached. We attribute this enhancement to a broad spectral function of antikaons at high 
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baryon density and to the 'subthreshold' energy of 1.8 A GeV considered. However, even 
in case of the off-shell propagation of baryons and mesons the experimental K~ spectra 
are underestimated and attractive antikaon potentials are still needed to achieve a proper 
description. 

At AGS energies (^ 11 A GeV) the particle production in the HSD approach essen- 
tially occurs via the excitation and decay of strings which can be viewed as continuum 
excitations of hadrons. Any spectral broadening of the 'continuum' thus is not likely to 
be seen in the asymptotic particle spectra of pions, kaons or antikaons especially since 
they are most abundantly produced far above the individual A^A^ or vrA^ thresholds. 

We finally point out that although our off-shell transport approach appears to be in 
a reasonable agreement with the differential experimental data at least for SIS energies 
(except for K~ when discarding antikaon potentials), the question of proper off-shell 
transition amplitudes in the collision terms remains an open problem that has to addressed 



in the near future. Some steps in this direction e.g. have been taken in Ref. [oC 
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Figure 1: upper part: Pio, Piz and Mj as a function of z{t) for a momentum-dependent 
imaginary potential with Cw = 0.6 GeV^ and Aw = 2.0 GeV (lower part). The vacuum 
width is chosen as Fy = 0.02 GeV^ and the initial separation in the mass parameter of 
the testparticles is AM = Fy/20. 
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Figure 2: upper part: PjO) Piz and M, as a function of z(i) for a momentum-dependent 
complex potential with Cy = -0.3 GeV^, Ay = 2.0 GeV, Cw = 0.6 GeV^ and Aw = 2.0 
GeV (lower part). For the vacuum width and the initial mass separation the same values 
are used as in Fig. 1. 
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as a function of the invariant energy y^ for different collisional width F of the antikaon 
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Figure 7: The number of baryon-baryon (BB) collisions as a function of the invariant 
energy ^/s for Au + Au at 1 A GeV and 6 = 1.5 fm. The solid line is obtained from 
including the off-shell propagation of baryons in the transport approach while the dashed 
line stands for the result in the on-shell limit. 
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Figure 9: (l.h.s.): The inclusive momentum spectra of positive pions at 6iab = 44 ± 4° for 
Au + Au at 1 A GeV in comparison to the experimental data from the KaoS Collaboration 
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(r.h.s.): The inclusive momentum spectra of positive kaons at 6iab = 44 ± 4° for Au + Au 
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circles |^l|, open circles [0). The dashed histograms show the results from the on-shell 
propagation while the solid histograms are obtained including the off-shell propagation 
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spectrum provide an estimate for the statistical error bars of the calculations. 
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Figure 11: The inclusive spectra of positive kaons for Ni + Ni at 1.8 A GeV including the 
off-shell propagation in the transport approach (solid histogram, full triangles) and in the 
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Figure 12: The inclusive spectra of K~ for Ni + A^i at 1.8 A GeV including the off-shell 
propagation in the transport approach (solid histogram, full triangles) and in the on-shell 
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Figure 13: The rapidity distributions of protons, vr"^, K'^ and K~ for central collisions 
of Au + Au at 11.3 A GeV in comparison to the data from Refs. [^, ^. The solid 
histograms are obtained including the off-shell propagation in the transport approach 
while the dashed histograms result from the on-shell limit. 
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Figure 14: The transverse mass spectra of 7r+ (l.h.s.), K~^ and K mesons (r.h.s.) for 



central collisions of Au + Au at 11.3 A GeV in comparison to the data from Refs. |p9 |. 
The solid histograms are obtained including the off-shell propagation in the transport 
approach while the dashed histograms result from the on-shell limit. 
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